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_ ■ Abstract 

C/} ; A conservation law and stability, recovering phenomena and characteristic patterns of a nonlin- 

ear dynamical system have been studied and applied to biological and ecological systems. In 
our previous study, we proposed a system of symmetric 2n-dimensional conserved nonlinear dif- 
ferential equations with external perturbations. In this paper, competitive systems described by 
2-dimensional nonlinear dynamical (ND) model with external perturbations are applied to popula- 
tion cycles and recovering phenomena of systems from microbes to mammals. 

The famous 10-year cycle of population density of Canadian lynx and snowshoe hare is nu- 
merically analyzed. We find that a nonlinear dynamical system with a conservation law is stable 
and generates a characteristic rhythm (cycle) of population density, which we call the standard 
rhythm of a nonlinear dynamical system. The stability and restoration phenomena are strongly 
related to a conservation law and balance of a system. The standard rhythm of population density 
is a manifestation of the survival of the fittest to the balance of a nonlinear dynamical system. 

1 Introduction 

The concept of stability is important in order to understand natural phenomena in biological and 
engineering systems. In our previous study, we studied the relation between a conservation law 
and stability of a 2n-dimensional competitive system that contains competitive interactions, self- 
interactions and mixing interactions. The system consists of 2n-dimensional nonlinear differential 
equations required from Noether's theorem [1]. The 2n-dimensional nonlinear ordinary differen- 
tial equations for a competitive system constructed to satisfy the conservation law have properties 
such as the addition law, which is empirically interpreted as recovery from injuries of skin and 
tissues in biological bodies. 

It has been shown by many researchers that a relatively simple set of interactions can ex- 
plain complex phenomena in biological systems ElO. For example, in 1952, Turing suggested 
chemical molecular mechanism called the reaction-diffusion system [4] which is defined as semi- 
linear parabolic partial differential equations. This reaction-diffusion system is well applied for 
explaining stripe patterns of the marine angelfish, Pomacanthus, and restoration phenomena in its 
stripe patterns from injuries was observed B21 |6l Q. Prigogine also proposed Brusselator model 
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with nonlinear ordinary differential equations to illustrate spatial oscillations and Turing patterns 
|[8l . It is also an interesting problem to investigate in ecological systems if a large complex sys- 
tem should be stable or not, and many researchers have discussed the criteria concerning the 
stability of a system for n dimensional ordinary differential equations and statistical framework 
[]9l [lOl dU [12l [131 EH H5). What would be a reason why a simple set of interactions can explain 
complex phenomena? We discussed a system of interactions generalizing Lotka-Volterra type non- 
linear competitive interactions and suggested that a conservation law could be a key to understand 
complex phenomena even in biological and ecological systems. 

We investigated the system of 2n-dimensional coupled first-order differential equations by us- 
ing Noether's theorem, which led to the following results, (i) The form of differential equations and 
coefficients of nonlinear interactions are strictly confined when the system has a conservation law 
which is constructed by interacting species of a particular experimental system, (ii) The conserved 
quantity of a system produces a Lyapunov function which is usually employed to study solutions 
of nonlinear differential equations. The conserved quantity is constructed by Noether's theorem, 
but the analysis of Lyapunov function would be used to check solutions to differential equations 
including those for non-conservative and dissipative systems. The system of differential equations 
with conservation law is different in this respect, (iii) A system of interactions could be analyzed 
as an assembly of a basic binary-coupled form (BCF). In other words, a complex interacting sys- 
tem can be decomposed into an assembly of binary-coupled systems. The BCF system is a simple 
basic set to explain complex phenomena defined by Noether's theorem, (iv) The BCF system with 
conservation law indicates an addition law which may be interpreted as the restoration or reha- 
bilitation phenomena; those are known in a large system of neural network or computer network 
when a small disordered device or a part of network system is replaced by a normal device. These 
properties could be applied to stability and restoration phenomena of biological systems, (v) The 
conservation law is also useful to check accuracy of numerical solutions to nonlinear differential 
equations. As summarized above, we discussed that the basic nonlinear system in BCF is stable. 
The binary-coupled system and addition law supported by a conservation law can lead to a large, 
stable complex system. This is an important conclusion in the conserved binary-coupled model. 
Because the BCF system has such several interesting properties, we will apply the model in order 
to study stability and interaction mechanism of a biological system. 

In this paper, we will explain the properties of solutions with a conservation law and applica- 
tions to biological systems. In Section [2l we extend the BCF model to simulate external pertur- 
bations numerically. There are various prey-predator type competitive models with perturbations 
[TT6l ITTl [T8l |T9l [20l [2T1 [22ll . however, most of them are with small, stochastic perturbations. The 
behaviors of conservation laws with external perturbations have been seldom considered. We will 
explicitly discuss properties of the conserved, stable, 2-variable nonlinear interacting system with 
external perturbations and the conservation law, its indications and possible applications to non- 
linear interacting system. We will show that 2-variable ND model has the properties of restoration 
and recovery from external perturbations. In Section [31 stability and population cycles of biolog- 
ical systems are examined in terms of a conservation law of the system. We will also examine 
specific examples of the Canadian lynx and snowshoe hare J23l |M |25l |26l W\ IM |29l and the 
question of population cycles, and food chain of microbes in the lake 11301 l3"TTl . Conclusions and 
summary of results are given in Section HI 
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2 The model of binary-coupled form (BCF) 



2.1 2n-ND system with perturbations 

We discussed BCF system and the conservation law of 2n-nonlinear dynamical (2n-ND) model in 
detail in the previous work [OQ. In this study, we include external perturbations in 2n-variable non- 
linear differential equations in order to examine characteristic behaviors of conserved nonlinear 
interacting systems. It should be noticed that the 2n-ND model is extended by adding external per- 
turbation terms which maintain a conservation law given by Noether's theorem. The odd variable 
terms for x,i (i = 1, . . . , 2ri) are 

n 

d-2k,2k-lX2k-l = ( a {2™+2fc} + &{2n 2 +2nk+2i-l})X2i-l + («{2n 2 +2m+2fc} + Ot{2n 2 +2nk+2i})X2i 

+ Ot{4n 2 +2ni+2k} x 2i-l%2i\ + OC{An i +2nk+i} x i x 2k-X + C2fc-1, 

3=1 

where k = 1, . . . , n. The even variable terms for Xi (i = 1, . . . , 2n) are 

n 

d-2k~l,2kX2k = | (Q{2m+2fc-l> + a {2nk+2i-l}) x 2i-l + (tt{2n 2 +2ni+2fc-l} + OL {2nk+2i}) X 2i 

(2-2) 

+ tt{4n 2 +2m+2fe-l} a; 2i-1^2i| + &{4n 2 +2nk+j}X jX2k + C2fc, 

where x = dx/dt, coefficients, d hj express d 2k ,2k-i = «2fc - a 2k ~i, d 2 k-i,2k = a 2k~i ~ a 2k - 
The linear coefficients and nonlinear coefficients a^, (i — 1, . . . , 8n 2 + 2n) are arbitrary constant 
values. The last terms c 2 fc-i, c 2 fc, (k = 1, . . .,n) of (12.11) and (12.21) are constant or piecewise 
continuous constant, which are interpreted as external perturbations (temperature, seasons and 
other temporal, external inputs). One should note that constant terms have dimension of velocity, 
so they are different from actual external perturbations which are considered to effectively express 
external perturbations. Because external perturbations (inputs) change population densities as x = 
dx/dt, we simulate numerically those effects with c 2 k-i, c 2 k as external inputs. The system has a 
conservation law derived from Noether's theorem which is proved in the paper JTJ: 

n 2n 

^ = ' < y^ J ^Ot{2ni+j}X2i-lXj + Ot{2n 2 +2ni+j}X2iX j + a{4n 2 +2ni+j}X2i-\X2iXj 

i=1 j=1 (2.3) 

n 

+ y^\c2jX2i-l + C2i-lX 2 i}. 
i=l 

Therefore, with the equations from (12.11) to (12.31) . we are able to consider the conserved nonlinear 
dynamical system with external perturbations by employing piecewise continuous constant terms, 

C2k-l, C2k- 

2.2 Properties of 2-variable ND model 

The equations of 2-variable ND model are produced by setting n — 1 (k — 1) in equations (|2.1I) to 
(12.31) . resulting in 

1 C\ 

x 1 = — {(a 4 + a 5 )xi + 2a 6 x 2 + 2a 8 XiX 2 + a 7 x\} + — , (2.4) 

"21 "21 
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1 C 2 

x 2 = —{2a 3 xi + (a 4 + a 5 )x 2 + 2a 7 xix 2 + a 8 xl} + — , (2.5) 

"12 "12 

and the 2-variable ND model has the following conservation law, 

= a 3 x x + {0L4 + a^)x\x 2 + o^a^ + «7X 1 X2 + ag^i^ + c 2 x i + c i^2- (2.6) 

The nonlinear interactions can generally represent, for example, Lotka-Volterra type prey- 
predator, competitive interactions, food-chain relations by adjusting nonlinear parameters ai, . . . , aj 
The piecewise continuous constants, c\ and c 2 are used as external perturbations in computer simu- 
lations, such as environmental conditions, food or hormones which increase or decrease interacting 
species in questions. The equations (12.41) ~ (12.61) form 2-variable BCF nonlinear differential equa- 
tions with a conservation law. 

By employing eqs. (|2.4I) ~ (12.61) . we will show: 

(1) solutions to the binary-coupled nonlinear equations maintain a characteristic (xi,x 2 ) phase- 
space of solutions and recovery from external perturbations. The external perturbations can nu- 
merically reproduce environmental conditions such as temperature and climate, environmental 
hormones and chemical substances which affect interacting species. The nonlinear binary-coupled 
model can be applied to examine responses of a system whether they are induced from internal 
interactions or external perturbations. 

(2) The binary-coupled nonlinear equations with conservation law exhibit stable phase- space solu- 
tions, which are interpreted as stability and recovery of population-change in a biological system. 
The properties of the binary-coupled nonlinear interactions will be shown explicitly in numerical 
simulations. 

(3) By employing the 2-variable binary-coupled model, it is possible to simulate cycles of maxima 
and minima in population-change, delays of periodic times of population cycles for competitive 
species [16]. Hence, cycles of population-change will be discussed in terms of the conservation 
law and nonlinear interactions. 



Figure |2.1(a)| shows the nonlinear interactions between species without external perturbations 
(ci = and c 2 = 0), whose coefficients of nonlinear equations are set as in Table 1 (Condition 1). 
In a view of the classical Lotka-Volterra competitive system, it can be interpreted as that x\ and 



x 2 represent prey and predator, respectively. Figure 2.1(b) is the phase-space given by solutions 



(xi,x 2 ). The solutions {xi,x 2 ) are periodic with respect to time, the maximum and minimum 



of (xx, x 2 ) appear with a time-delay. Figure [27l(c)| shows the numerical value of the conserved 



function ^ defined by (12.61) . which is constant with respect to time. 



The solutions (xi, x 2 ) in Figure 2.1(a) show explicitly a time-delay of the peak for interact- 



ing species. The timings of peak and delayed peak are determined by nonlinear interactions and 



strength of coupling constant. The solutions [Xi, x 2 ) in Figure [271 (b)| show phase-space solutions, 
which are stable in the meaning that the conserved quantity \& is maintained constant and phase- 
space solutions are in the same trajectory for all time. The unit of time should be considered to 
adjust to the time scale of a system in consideration, because biological unit times are generally 
different from microbes to mammals. 



Table 1: The list of nonlinear coefficients. 
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(a) 2-variable ND solutions. Solid and dashed lines represent x\ 
(prey) and X2 (predator), respectively. One should note that the 
unit of time should be defined with respect to a system in consid- 
eration. 
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(b) Phase-space of 2-variable ND solutions. (c) Conservation law ^ of 2-variable ND. It is constant with re- 

spect to time. 



Figure 2.1: A 2-variable ND solution and Conservation law ^. 



The phase-space diagram [2. l(b)| and the straight line of Figure [271 (c)| show that the solution is 
exact and stable [1]. The three figures exhibit important properties of solutions to the system of 
prey-predator type of competitive nonlinear interactions. 

One of the important properties shown by the stable, conserved nonlinear system is that the 
interacting species repeat the rhythm of maxima and minima of the population. The periods of 
the rhythm are the result of complicated nonlinear interactions, but the system keeps the constant 
quantity ^ with respect to time. The interesting applications of the BCF model are shown by em- 
ploying in the paper "My sis in the Okanagan Lake food web ' QUI . Canadian Lynx and snowshow 
hare [24), which will be explained in Section 3. 

2.3 Recovering and restoration from perturbations 

In order to investigate responses of a system to external perturbations, we introduce piecewise 
continuous constants, c\ and ci, by using ^-functions such that 

Q = fi{6{t - t start ) - 6{t - t end )}, (i = 1,2), (2.7) 

where 6(t — £') represents a step function: 



e{t-t') = {' (2.8) 




and coefficients f if (i = 1,2) are positive or negative constants to express strength of external 
perturbations. The constants are adjusted to produce reasonable maxima and minima in numerical 
simulations. 
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(a) 2-variable ND solutions with a negative perturbation on prey 
si. The perturbation is introduced from t = 700 to t = 1200 
which is represented as gray background. 
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(b) The (xi,X2) phase-space transition with the negative per- 
turbation as in (a). Solid line (St. 1) is initial state, and St. 2 is 
the recovered state after the end of perturbation. Dashed line is 
phase-space during Sp. 1 - Ep. 1. 
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(c) Conservation law * of 2-variable ND with one perturbation. 
^ changed * ~ 60000 to * ~ 30000 by introducing perturba- 
tion. It recovers after Ep. 1 . 



Figure 2.2: An external perturbation and a recovery. 



Figures |2.2(a)[ |2.2(b)| and |2.2(c)| show the reaction and recovery of the nonlinear interacting 
system from an external perturbation. One of the typical recovery of a system from a perturbed 
state is shown. In Figure [2~2(a)[ an external perturbation starts at t = 700 (Sp. 1), and the coefficient 
fi equals to —1260.0 and / 2 equals to zero in this example. The black arrow is the starting point 
of perturbation, and the gray arrow is the end of perturbation in Figures |2.2(a)| and |2.2(c)[ The 
nonlinear coefficients are listed in Table 1 (Condition 1). The solutions (xi, x 2 ) are deformed by 



the perturbation (Figure [2.2(a)| and |2.2(b)[ ). However, the system does not disintegrate but finds a 
new stable phase-space close to the original phase-space and maintains a new conserved relation. 
The perturbation ends at t = 1200 (Ep.l), and the system recovers the original state (x\, x 2 ). 

The timing of negative perturbation which reduces the population number x\ or x 2 produces 
different results. When a negative perturbation is exerted in the increasing phase of xi or x 2 , the 
system will find a new conserved stable solution near the original solution, but when a negative 
strong perturbation is exerted before x\ or x 2 gets to its minimum, the system may collapse: the 
system exhibits no solutions (x±, x 2 ), which would be interpreted as disintegration or extinction in 
biological systems. 
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(a) 2- variable ND solutions with a critical negative perturbation 
on prey, x\. Solutions converge to zero after the perturbation. 
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(b) Conservation law \l> with a critical negative perturbation. ^ 
converges to zero after the critical perturbation. 



Figure 2.3: Critical negative perturbation and extinction. 



The conserved nonlinear system naturally exhibits maxima and minima without external per- 
turbations, and so we call these maxima and minima as endogenous maximum and minimum. It is 
needed to distinguish them from enhanced maxima and minima by external perturbations. 

In Figure I2.3L the response of a strong negative perturbation to prey after the peak of endoge- 
nous maximum is shown. The values of coefficients are listed in Table 1 (Condition 1). The 
starting point of this perturbation is at t = 800 and the end point of the perturbation is at t = 950. 
The negative constant of perturbation is f% = —3175.3879. The prey, X\, rapidly declines with 
negative perturbation, and (x\, x 2 ) converges to zero for t > 1000. These computer simulations 
may be compatible with known empirical results, for example, in pest control. A pest control 
is not so effective if it is performed in the season when harmful insects are in peak and active, 
because species are energetic enough to find a new stable life to live. It is effective when a pest 
control is performed in the season when harmful insects are not so active or in a declining state 
after endogenous maximum. 

In the nonlinear interacting system, positive perturbations which will increase x\ or x 2 do not 
always mean a positive effect on stability of the system. There is a limit to the value of a positive 
perturbation, because an increase of x\ leads to a decrease of x 2 in a stable system (xi, x 2 ), which 
indicates that a system has internally allowed maximum and minimum populations. 

Figures I2~4l shows the behaviors of (x 1 , x 2 ) at normal and critical values of positive perturba- 
tions, ci, for x\. The values of coefficients are listed in Table 1 (Condition 2). Figures |2.4(a)| and 



2.4(b)| show m at the normal positive perturbation which increases interacting species will increase 



the peak of (xi,x 2 ) populations. However, at certain critical values of coupling constants, the 
prey-predator interaction cannot keep and support the rhythm of maxima and minima, and the sys- 
tem diverges. Figures [2 .4(c) | and |2 .4(d)| show that the system cannot maintain a stable, interacting 
system when the positive perturbation surpasses the critical value (ci = 1599.924999 in the current 
simulation). The unstable solutions branch out at t ~ 1100 when the value of perturbation changes 
from f x = 1160.0 to f\ = 1599.924999. 

Hence, in a conserved stable system, species seem to strictly control each other by seeking a 
new stable solution so that they can survive together. The competitive interacting system such as 
the conserved prey-predator relations may be considered to be a cooperative system for species to 
survive. It should be noted that if a dynamical prey-predator system is active, the rhythms of max- 
ima and minima are clearly repeated, which is known in real prey-predator systems. However, if 
an external perturbation (exogenous interaction) exceeds a certain critical value of the competitive 
system, the rhythms of maxima and minima will disappear first and then after a time, the system 
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(a) 2-variable ND solutions with a positive perturbation. The 
perturbation starts at t = 500 and ends at t = 900 represented 
as gray background. The amplitudes of x\ and X2 become larger 
than before. 
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(b) The Conserved quantity with a positive perturbation. It re- 
covers after perturbation but finds another equilibrium state. 
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(c) 2-variable ND solutions with a critical perturbation. x\ and (d) Conservation law of 2-variable ND with a critical perturbation. 

X2 converge to zero after a critical perturbation. The ^ converges to zero after a critical perturbation. 



Figure 2.4: Critical perturbations and divergence of solutions. 



will diverge (disintegrate). Therefore, the rhythm of wild-life indicates that the dynamical inter- 
actions between species are active and stable. When the rhythm of change disappear or doesn't 
come back, it may indicate that related species are in danger of extinction. The rhythm is impor- 
tant to examine if the wild life is normal and active, or harmed by human activities and external 
perturbations. 

On the other hand, by adding another perturbation, we can show that it is possible to save 
species from extinction. Figure [275(a)| is a result of a positive perturbation to save species (xi, x 2 ) 
in a danger of extinction in Figure [273(a)| We exerted a positive perturbation after Sp. 1 - Ep. 1 in 
Figure [231 The positive perturbations start at t = 1000 (Sp. 2) and end at t = 1300 (Ep. 2), the 
strengths of c\ and C2 are fi = 200, f 2 = —1000. Figure [2~5 (a) | show s that species are in danger of 
extinction, however, if positive external perturbations are properly inserted, the system will come 
back to life again. 



3 Conservation law and population cycles 
3.1 The food- web of Microbes in Okanagan Lake 

One of interesting data of the ecological interactions is the interaction described in l Mysis and 
the other zooplanktons in the Okanagan Lake' OOl . Although the food- web in Okanagan Lake is 
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(a) 2-variable ND solutions with perturbations to avoid con- (b) Conservation law ^ with two perturbations in (a). ^ recovers 

verging to zero after a critical perturbation. x\ and X2 come from the perturbation after Sp. 2 - Ep. 2. 
back to life after the second perturbations. 



Figure 2.5: The critical behavior and restoration. 



not clarified definitely, mysis introduction to lakes is known as an effective method to enhance 
ecological interactions and its strengths among microbes and other creatures so as to increase 
fisheries productions. 

The time-series of dominant crustacean zooplankton densities in Okanagan lake has been mea- 
sured monthly and suggested that mysis and zooplankton populations are synchronous and charac- 
terized by the cycle of the peak and bottom population densities. The cycles of population densities 
are primarily due to cycles of season and climate and then to mutual interaction of microbes. The 
analysis of microbes suggests that the density-dependent and delayed population regulation of 
microbes is evident. In addition to the seasonal factors, the regular cycles and the delayed peak 
and bottom populations densities of microbes are the results of strong nonlinear interactions of 
species. We numerically examined changes of population densities of microbes by employing the 
2-variable conserved ND model. 

The current conserved nonlinear model shows that the interacting species designate a standard 
rhythm of the peak and bottom population densities. There are some fluctuations at the peak 
and bottom densities, but they show the stable dynamic life as demonstrated in Figure [2761 (a) ~ 
(c). Although, normal peak and bottom densities can be readily explained by adjusting coupling 
strength of model's internal interactions, a sudden change of maxima which is often encountered 
in a biological data cannot be easily simulated by only adjusting internal coupling constants in the 
2-variable nonlinear interacting model. 

In Figure I2.6L several perturbations are exerted on the interacting 2-variable system. The first 
external perturbation starts at t = 500 (Sp. 1) and ends at t = 1000 (Ep. 1). The strength of 
perturbations in Sp. 1-Ep. 1 are fx = —800, f 2 = —100. The second external perturbation starts 
at t — 1400 (Sp. 2) and ends at t — 1900 (Ep. 2). The strengths of perturbations in Sp. 2-Ep. 2 
are fx = —50, f 2 = —120. The third external perturbation starts at t = 2200 (Sp. 3) and ends at 
t = 2600 (Ep. 3). The strength of perturbations in Sp. 3-Ep. 3 is set as fx = —500, f 2 = —50. The 
lines (xx, x 2 ) may represent for instance, the prey -predator interactions, species of food-chain, and 
species interacting with its environmental factors (temperature or some environmental hormones). 
Black arrows are starting point of perturbations, and gray arrows are the end of perturbations; 
parameters are listed in Table 1 (Condition 1). The time period is within t = 4000, initial values 
are Xx = 500, x 2 = 300. 

The significant properties of the stable nonlinear conserved system are that if external pertur- 
bations are not large enough to disintegrate the system, the system will find a stable conserved 
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(a) 2-variable ND solutions with three external perturbations. The 
rhythm of x\ and X2 recovers from several perturbations. Gray 
backgrounds represent periods of perturbations. 




(b) Phase-space transitions of x\ and X2- Dashed lines rep- 
resent solutions, xi and X2, during perturbations. Solid line 
represent solutions without perturbations. 
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(c) Conservation law with three perturbations. It recovers from 
three perturbations. * ~ 60000 in the St. 1 and St. 2. 



Figure 2.6: Several external perturbations and recoveries. 



solution near the original system and continue a stable cycle (maxima and minima). It is clearly 
seen from (xi, X2) -phase space solutions in Figure |2~6(b)| The system recovers from several exter- 
nal perturbations. 

The numerical analysis can be applied to examine the change of population densities of mi- 
crobes. For example, the time-series data of dominant crustacean zooplankton densities in the 
Figure 2 of the paper 'My sis in the Okanagan Lake food- web ', show that the sudden maxima of 
dominant zooplankton densities are seen in the period '99 ~ '02. The sudden increase of the peak 
is readily adjusted when an external perturbation is assumed in the simulation, however, it is not 
reproduced by adjusting internal coupling constants in the 2-variable nonlinear model. Hence, it is 
concluded in the 2-variable model that there would have been certain positive external perturbation 
to the system of microbes in Okangan Lake during '98 ~ '01 considering a time-delay of external 
perturbations. 

It is interesting to check what kind of external or internal perturbations is affecting the peak 
of population density during the period '98 ~ '01. If there are no explicit changes in external 
or internal factors during the period, a sudden increase of the peak could be a result of more 
complex internal interactions. For example, the rhythm of the peak and bottom population densities 
should be explained by 4-variable or 6-variable nonlinear interactions of microbes. The unusual 
rhythm indicates how exogenous (environmental) and endogenous (internal interactions) variables 
are affecting the dynamics of each component and environmental nature related to the species. The 
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analysis of nonlinear model suggests that the sudden peak and bottom densities have important 
information on the dynamics of the system of species and environment. Hence, it is important to 
understand the standard rhythm of the peak and bottom population densities in order to distinguish 
them from unusual maxima and minima. 

One should be careful that a positive perturbation on one of interacting species not only en- 
hances the peak of maxima but also decreases minima in the rhythm of species. It is often true 
that the effect of enhancement is usually emphasized without taking care after negative effects. 
Hence, the enhancement of the number of population of a specific species may be harmful to other 
species in the food-web and consequently it endangers itself. Our analyses in Figure [2~4l and [231 
show that if we carefully control the increase or decrease of the population of certain species after 
introduction of a positive effect, we can keep normal and stable dynamics of species suitable for 
the environment. For this purpose, it is essential to explicitly understand the standard rhythm from 
real observed data. 

3.2 Population regulation in Canadian lynx and snowshoe hare 

It is difficult to identify population regulation mechanisms about prey-predator patterns of large 
mammals because the large mammal's life span is relatively long compared with microbes. The 
prey-predator cycle such as wolves and caribous takes some decades of years to observe, their in- 
teracting relation and behaviors have been recently revealed with modern technology (GPS-colored 
animals) ||23~1 . However, the food- web configuration between snowshoe hare and Canadian lynx is 
well-known prey-predator type phenomena, and a ten-year cycle of Canadian lynx was examined 
from the data of Canada lynx fur-trades return of the Northern Department of the Hudson's Bay 
Company (the data are from C. Elton and M. Nicholson [24]). 

The Canadian lynx and snowshoe hare have a synchronous ten-year cycle in population num- 
bers Il23l 1251 . The fundamental mechanisms for these cycles are maintained by the important 
factors such as nutrient, predation and social interactions [26J. In addition to the important factors, 
the nonlinear model with conservation law suggests that species of a system consequently find a 
strategy or a mechanism to survive for long-time periods. In other words, the cycle of population 
density is a manifestation of the strategy or mechanism to survive, which is suggested by stability 
of phase-space solutions determined by conservation law of a system. 

The nonlinear interactions with conservation law show a standard rhythm and stability from 
external perturbations as shown in Figure 12.61 The feeding and nutrient experiments in 11261 are 
considered as external perturbations to the system. As shown in Figure I2.6L the perturbations 
cause certain effects on the system, but the system will find a rhythm to maintain the dynamics of 
species, which is not so different from the original standard rhythm. Our numerical results agree 
with conclusions derived from feeding experiments and nutrient- addition experiments. Therefore, 
we propose that the properties of the system which has a conservation law should be a key to 
understand the unanswered question: why do these cycles exist?. 

The results of computer simulations show that the timing of perturbation leads to different 
results. This is also confirmed by the feeding experiment of snowshoe hare: " ... during the 



Table 2: The list of nonlinear coefficients of simulation in Figure [3TT1 
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a 3 
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1.0 


635.0 


-0.35 


700.5 


300.5 


-0.35 


-0.0068 


-0.016 
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(a) The 2- variable ND simulation of Canadian lynx population. 
The solid line represents Canadian lynx population | 24), and 
the dashed line represents a theoretical solution of 2-variable 
ND with several perturbations. 
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(b) The estimated population of Canadian lynx and snowshoe 
hare. The dashed line represents Canadian lynx population sim- 
ulated by 2-variable ND model with perturbations, and the solid 
line represents approximate population of snowshoe hare. 



80 100 



Time 

(c) Transition of conservation law * with respect to time. Several 
perturbations are introduced. 



Figure 3.1: Simulation of Canadian lynx and snowshoe hare. 



peak of the cycle in 1989 and 1990 had no impact on reproductive output ... however, during 
the decline phase in 1991 and 1992, the predator exposure plus food treatment caused a dramatic 
increase in reproductive output ..." Eoll . This fact can be examined in our model calculations. The 
perturbation in the peak phase does not cause large effects on standard rhythm, but negative and 
positive perturbations during a decreasing or increasing phase induce dramatic effects. 

The cycle of standard rhythm for Canadian lynx and snowshoe hare indicates that the stable 
dynamical system of lynx and hare functions normally and environmental nature is conserved 
in reasonable conditions. However, as we have shown in Figure |2.4(c)| and |2.4(d)[ if a strong 
negative perturbation is applied persistently for a long period, the system would fall into a danger 
of extinction. The important results of our simulation tell that before a system gets in danger of 
extinction, the standard rhythm of the system will tend to become ambiguous or disappear. Hence, 
if we carefully observe the standard rhythm of a specific system of species, we could help the 
dynamical system save and preserve related natural environment. 

In FigureED we simulated the Canadian lynx data of the Hudson's Bay Company from 1821 to 
1910, which is approximately thought as the lynx-population density. The interpolated Elton's data 
was downloaded from lhttp : //www, atomosycl. net/spip . php?action=dw2_out &icl=4"2| 



The solid-line in Figure [371 (a) \ is lynx-population data and the dashed-line is the results of our nu- 
merical simulation using 2-variable nonlinear interactions between lynx and snowshoe hare (Figure 



3.1(b)[ ). The conserved binary-coupled model tells that there should have been some external per- 
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Table 3: The list of external perturbations in Figure [3TTT The periods of positive and negative 
perturbations to numerically simulate Canadian lynx population. Note that the values of fi have 
the meaning of velocity (number/time). 



Perturbation 


Time 


Strength of /i 


Sp. 1 - Ep. 1 


<t< 11 


-11000000 


Sp. 2 - Ep. 2 


20 < t < 30 


-4000000 


Sp. 3 - Ep. 3 


30 < t < 35 


-12000000 


Sp. 4 - Ep. 4 


35 < t < 41 


-8000000 


Sp. 5 - Ep. 5 


41 < t < 46 


4200000 


Sp. 6 - Ep. 6 


50 < t < 60 


-9000000 


Sp. 7 - Ep. 7 


60 < t < 69 


1140000 


Sp. 8 - Ep. 8 


70 < t < 80 


-11000000 


Sp. 9 - Ep. 9 


87 < t < 100 


-14500000 


Sp. 10 - Ep. 10 


100 < t < 113 


21500000 



turbations, although we cannot make sure at the present what kinds of external perturbations were 
exerted. The actual population density of snowshoe hare is not known, and so we assumed a rea- 
sonable population density and several external perturbations for numerical simulations in order to 
fit the lynx population data (see, Table 2 and Table 3). 

The snowshoe hare gets several positive and negative perturbations, but the overall rhythms 
of lynx and hare are not altered. As suggested by in Figure |2.6(b)[ the phase-space of lynx and 
snowshoe hare is stable against several external perturbations. This is also compatible with the em- 
pirical fact that the ten-year cycle in snowshoe hare is resilient to a variety of natural disturbances 
from forest fires to short-term climatic fluctuations. However, as shown in our model calculation 
in Figure [274(c) [ a long-term (more than ten years) negative perturbations and a vast environmental 
change that humans could cause would definitely endanger the standard rhythm of snowshoe hare, 
lynx and related species. 



4 Conclusions 

In this paper, we examined characteristic properties of several ecological systems based on con- 
served nonlinear interactions which include generalized Lotka-Volterra type prey-predator, com- 
petitive interactions. In Section |2~TI we extended our 2ri-variable ND model by including external 
perturbations in order to apply the model to more realistic biological phenomena and to study 
responses of a biological system to external perturbations. 

We simulated external positive and negative perturbations by employing piecewise constant 
terms in our nonlinear equations. As it is discussed in the analysis, the results of simplified per- 
turbations agreed with the experiments and empirical data reasonably well. The numerical simula- 
tions showed the existence of the standard rhythm which is characteristic to a nonlinear conserved 
system. It is essential to understand standard rhythm by observing and taking data of a system so 
that we can distinguish unusual maxima and minima from standard rhythm. This gives a possibility 
to examine signatures that distinguish internal effects from external ones. 

The ten-year cycle of lynx and hare is a very interesting biological phenomena. Though a 
cycle of a biological system should be a phenomenon composed of complex and multi-biological 
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interactions, the 2-variable BCF analysis has revealed the interesting results on properties of the 
biological phenomena. The ten-year cycle of lynx and hare is stable and resilient to external 
perturbations, which is reproduced in our model calculations. The system with conservation law 
shows stable cycles and recovering phenomena, which are displayed numerically in phase-space 
solutions. The stability and conservation law are constructed at least by binary-coupled species in 
biological and ecological systems, and they are maintained in a more complicated multi-coupled 
system, as we proved in a general form (TJ. 

The coupling constants of interacting species expressed in nonlinear differential equations are 
considered to have been determined in a long time by complicated environmental and internal fac- 
tors of a specific system, such as the landforms, seasons, climate and temperature. Once members 
and structures of dynamical systems were constructed, appropriate dynamical systems would be 
maintained for long-time periods with internal factors such as nutrient, predation and social inter- 
actions. The predation and social interactions are expressed as complicated nonlinear relations in 
mathematical terms. This may be explained by the fact that members of a system have a well- 
conserved rhythm respectively and these rhythms also have a well-determined slight delay to each 
other, which indicates that certain nonlinear interactions among members exist. 

The important factors (nutrient, predation and social interactions) are needed for all species to 
survive in nature, but they easily change by natural conditions. In addition, an unusual increase of 
population numbers of a species would endanger the survival of a species itself as well as other 
species (see the numerical simulations in Figure 12.41) . The important property of the nonlinear 
model with conservation law is that the binary-coupled system can have the persistent stability and 
recovering strength to external perturbations. As a predator needs a prey for its food, a prey needs 
a predator for the conservation of their own species. The conservation law and rhythm of species 
are considered to be constructed by species and natural conditions in a system for a long time, and 
hence, the cycle (rhythm) of species would be interpreted as a manifestation of the survival of the 
fittest to the balance a biological system. 

We conclude that stability and conservation law are constructed by species in mutual depen- 
dency or cooperation to survive for long-time periods in severe nature. The standard rhythm should 
be regarded as the result of strategy for species to live in nature. Whatever roles they have to play, 
the species that can fit and balance with other creatures can survive in nature. A strong predator 
cannot even survive if it ignores the law of the standard rhythm and conservation law of a system 
constructed by other members and the environment. We hope that this study will help understand 
both activities of animals and humans in natural life. 
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